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ABSTRACT 

We re-examine the assembly history of field spheroidals as a potentially powerful dis- 
criminant of galaxy formation models. Whereas monolithic collapse and hierarchical, 
merger-driven, models suggest radically different histories for these galaxies, neither 
the theoretical predictions nor the observational data for field galaxies have been suf- 
ficiently reliable for precise conclusions to be drawn. A major difficulty in interpreting 
the observations, reviewed here, concerns the taxonomic definition of spheroidals in 
merger-based models. Using quantitative measures of recent star formation activity 
drawn from the internal properties of a sample of distant field galaxies in the Hub- 
ble Deep Fields, we undertake a new analysis to assess the continuous formation 
of spheroidal galaxies. Whereas abundances and redshift distributions of modelled 
spheroidals are fairly insensitive to their formation path, we demonstrate that the 
distribution and amount of blue light arising from recent mergers provides a more 
sensitive approach. With the limited resolved data currently available, the rate of 
mass assembly implied by the observed colour inhomogeneities is compared to that 
expected in popular A-dominated cold dark matter models of structure formation. 
These models produce as many highly inhomogcneous spheroidals as observed, but 
underpredict the proportion of homogeneous, passive objects. We conclude that colour 
inhomogeneities, particularly when combined with spectroscopic diagnostics for large, 
representative samples of field spheroidals, will be a more valuable test of their physi- 
cal assembly history than basic source counts and redshift distributions. Securing such 
data should be a high priority for the Advanced Camera for Surveys on Hubble Space 
Telescope. 



1 INTRODUCTION 

The star formation and assembly history of field 
'spheroidals' (defined here to include both ellipticals and 
SOs) has been discussed as one of the most powerful tests of 
galaxy formation capable of distinguishing between very dif- 
ferent structure formation possibilities (Kauffmann, White 
& Guiderdoni 1993; Baugh, Cole & Frenk 1996b; Kauff- 
mann & Chariot 1998a). 

Traditionally spheroidals were considered to have 
evolved as isolated systems following a 'monolithic' collapse 
at high redshift (Eggen, Lynden-Bell & Sandage 1962). Ev- 
idence in support of this conjecture comes from many quar- 
ters, including their high central stellar density and slow 
rotation, consistent with non-dissipative collapse at early 
epochs (Efstathiou & Silk 1983), and their remarkable ho- 
mogeneity at various look-back times in colour-magnitude 
(Sandage & Visvanathan 1978; Bower et al. 1990; Ellis et 
al. 1997; Stanford, Eisenhardt & Dickinson 1998) and fun- 



damental plane studies (van Dokkum et al. 2001; Treu et 
al. 2001) of distant clusters. 

This picture contrasts with that naively expected in 
models where galaxies assemble hierarchically at a rate gov- 
erned by the merging of cold dark matter halos and dy- 
namical friction timescales (Kauffmann, White & Guider- 
doni 1993; Baugh, Cole & Frenk 1996b; Cole et al. 2000). 
Distinguishing between both hypotheses is central to veri- 
fying the cold dark matter picture but progress is hindered 
by the need to distinguish between star formation histories 
and those of mass assembly. 

Most of the early work concentrated on ellipticals in rich 
clusters where there are significant difficulties of interpreta- 
tion. The early epoch of star formation consistent with the 
colour-magnitude and fundamental plane studies may still 
be consistent with hierarchical models since clusters form 
from the rare, high density peaks in the primordial fluc- 
tuations within which we can expect accelerated evolution. 
Moreover, comparisons of galaxies in individual clusters se- 
lected at various redshifts may be fundamentally affected by 
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the way in which rich clusters are themselves selected. Con- 
veniences of observation (e.g. selecting clusters of a given 
richness, with a relaxed distribution or of a given X-ray tem- 
perature) may bias us towards the selection of more evolved 
systems at a given epoch (Kauffmann 1996), thereby giv- 
ing a misleading impression of a uniformly old population. 
Most fundamentally, if galaxies within clusters are them- 
selves continuously evolving or arriving, it will be difficult 
to conduct all-inclusive inventories of a given morphological 
type as a function of look-back time. 

For the above reasons, attention has now turned to- 
wards providing observational constraints on the evolution 
of field galaxies. The term 'field' here includes systems in 
clusters, of course, and implies an attempt to be all-inclusive 
in the study of the population. The simplest test for evolu- 
tion in field spheroidals (and the one which has received 
most attention) is concerned with establishing their number 
density at high redshift, either using HST morphologies to 
isolate spheroidals (Zepf 1997; Schade et al. 1999; Frances- 
chini et al. 1998; Im et al. 2001), optical-infrared colours 
to locate passively-evolving examples whose stars formed at 
high redshift (Barger et al. 1999; Daddi et al. 2000; Mc- 
Carthy et al. 2001) or some combination of both (Menan- 
teau et al. 1999). 

We will show in this article that tests which rely on 
abundance comparisons at various redshifts or apparent 
magnitude limits, are currently inconclusive and, more fun- 
damentally, are poorly-suited to discriminating between the 
physical details of the various assembly hypotheses. In terms 
of the observations, the currently available spectroscopic 
datasets insufficiently sample the redshift domain, 1< z <2, 
where evolution is expected in popular low density models 
(Daddi et al. 2000; McCarthy et al. 2001; Firth et al. 2001). 
Comparisons based on counts alone are reliant on a good 
knowledge of the shape of the local luminosity function and 
its evolution. Without redshift data, there is considerable 
ambiguity in the implications (Glazebrook et al. 1995). In 
the case of samples defined morphologically from HST data, 
there is further confusion as to the merits of including com- 
pact objects which may possibly be unrelated star-forming 
galaxies whose abundance seems to increase significantly at 
faint magnitudes. Finally, selections based on red colours 
alone, even those utilising the new generation of infrared sur- 
veys, permit inclusion of unrelated dusty later-type galax- 
ies. Colour and morphology are known to poorly correlate 
at low to intermediate redshifts (Schade et al. 1999) and 
beyond z ~1 where dusty sources are more likely, the po- 
tential for confusion can only increase. The main drawback 
of abundance comparisons is that only subtle changes are 
expected at the depths currently being probed in low den- 
sity universes, thus the significance of claims for or against 
the hierarchical and monolithic (or number-conserving pure 
luminosity evolution) models remains marginal. 

Significant uncertainties are also present in the theo- 
retical models used to interpret the observations. In semi- 
analytic models the build up of the stellar mass of galaxies 
is determined using simple arguments for cooling and star 
formation rates in a merging hierarchy of dark matter ha- 
los. Even granting this very simple prescription, there is a 
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Figure 1. Theoretical uncertainties in the evolutionary history of 
massive spheroidals: the solid line and circles show the comoving 
density of all galaxies with stellar mass greater than IO^'^Mq as 
a function of redshift in the ACDM models of Cole et al. (2000) 
and Kauffmann et al. (1999) respectively. 

surprising uncertainty concerning the expected evolution in 
the number density of massive systems. For example, we 
can consider the comoving abundance of objects selected by 
stellar mass as a function of redshift. In Fig. ^ the solid line 
and circles show the comoving abundance of galaxies with 
stellar masses greater than 10^^ M0 in a ACDM cosmology 
as predicted by Cole et al. (2000) and Kauffmann et al. 

(1999) respectively. (Note that while the basic cosmological 
parameters, Qo and Ao were the same for these two calcu- 
lations others, including fib, were different. The Cole et al. 

(2000) results are taken from the models used in this paper, 
while those of Kauffmann et al. (1999) were computed from 
their publically available galaxy catalogues.) Although the 
two models are in reasonable agreement at 2: = only rather 
modest evolution (about a factor three decline) is predicted 
to 2 ~1 by Kauffmann et al. (1999), whereas much stronger 
evolution is expected according to the model of Cole et al. 
(2000). 

This discrepancy in the rate of evolution predicted by 
two similar models is very worrying, although it should be 
remembered that the galaxy stellar mass function is steep 
in this region. Hence even minor differences in the rate at 
which galaxy masses evolve in the two models could lead to 
large differences in abundance. For example, the factor of 10 
difference ai z — 1 might be resolved by considering galax- 
ies only a factor of 2 less massive in the Cole et al. model. 
Nevertheless, it is evidently crucial to understand these dif- 
ferences before strong tests of the hierarchical hypothesis 
can be made. 

A recent analysis of spheroidals in both Hubble Deep 
Fields (HDF) by Menanteau, Abraham & Ellis (2001) offers 
an interesting opportunity to bypass many of the difficulties 
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discussed above. Resolved V — I colour data for the HDF 
spheroidals reveals the presence of blue cores in a significant 
subset and preliminary spectroscopy (c.f. Ellis 2000) con- 
firms that recent star formation is the likely cause. Together 
with a careful appraisal of the way in which spheroidal galax- 
ies are defined in semi-analytic models, further analysis of 
this dataset forms the basis of the present paper. Instead 
of making absolute count comparisons, we use the distribu- 
tion of colour inhomogeneities to determine the fraction of 
spheroids still assembling. We compare this with hierarchi- 
cal predictions taking into account the selection technique 
and uncertainties arising from the taxonomical definition of 
a recognisable spheroidal galaxy. This type of comparison, 
in which the information of interest is encoded in a complex 
way in the data and important selection criteria exist, is one 
of the key areas where semi-analytic models excel. Most sig- 
nificantly, we demonstrate that such an analysis, based on 
post-merger relaxed systems with residual star formation, 
is a more direct test of the pertinent assembly details than 
traditional techniques based on counts or luminosities alone. 

A plan of the paper follows. In ^ we discuss the var- 
ious ways in which spheroidals can be defined in hierar- 
chical models, introducing a definition that incorporates a 
timescale for relaxation to a spheroidal form. We test pre- 
dictions utilising these definitions against published mor- 
phological number counts (Glazebrook et al. 1995; Driver 
et al. 1995a, 1995b; Franceschini et al. 1998) and redshift 
distributions (Brinchmann 1999; Menanteau et al. 1999). 
In §^ we revisit the observational analysis of Menanteau et 
al., and compare theory with observation in the light of both 
observational and theoretical uncertainties. We summarise 
our conclusions in §0. 



2 SPHEROIDAL GALAXIES IN 
HIERARCHICAL MODELS 

2.1 Semi- Analytic Modelling and Morphological 
Evolution 

First we discuss the basic processes which govern the mass 
assembly and star formation histories of spheroidals in hier- 
archical models. Specifically, we will utilise the galax;y for- 
mation model described by Cole et al. (2000) to which the 
reader is referred for a full description. Here we summarise 
briefiy the salient features. A small number of model pa- 
rameters are crucial in determining the properties of the 
spheroids in which we are interested as described below. 
Cole et al. describe how these model parameters may be 
constrained by a set of z = observational data. We retain 
the same parameter values which Cole et al. found produced 
the best agreement with the local data. As such, the model 
results presented in this work are true predictions. In §^ we 
will consider the effects of changing some of these param- 
eters (and also in altering the basic mechanisms by which 
spheroidal galaxies are made) as a way to understand the 
model results and to test their robustness and also to ascer- 
tain the prospects of our approach for directly constraining 
the formation history of spheroidals. It must be kept in mind 
however, that once a parameter value is changed the model 



is unlikely to remain a good fit to the local data and so 
any such results cannot be considered to be predictions of a 
realistic hierarchical model. 

In this model, galaxies form continuously from cooling 
gas inside a dark matter halo. Gas settles into a disk configu- 
ration within which stars form quiescently. Spheroidal galax- 
ies are produced by a major merger of two pre-existing galax- 
ies of comparable mass (note that the pre-existing galax- 
ies may be disks, spheroids or a mixture of both, and that 
disks may contain both stars and gas), which occurs due 
to dynamical friction (the timescale of which is determined 
by a parameter of the model, fdi, which is set to 1 unless 
otherwise stated) dissipating energy from the galaxies' or- 
bits. Such a merger disrupts the galaxies leaving a spheroidal 
remnant, and also triggers any cold gas present to undergo 
a burst of star formation with an exponentially decreasing 
star formation rate. The spheroidal so formed may accrete 
gas through cooling at a later time and thereby grow a new 
disk. Model galaxies can thus migrate either way along the 
Hubble sequence, and may change their morphologies sev- 
eral times during their lifetimes. 

Only a major merger (defined as one in which Mi /M2 > 
/eiiip, where Mi and M2 are the masses of the merging com- 
ponents with Ml < M2 and /eiup is a parameter of the 
model) triggers a burst of star formation. The timescale of 
the burst is Tburst = el'.burst'^dyn, where Tdyn is the dynami- 
cal time of the newly formed spheroid and e*, burst is a free 
parameter. 

Cole et al. (2000) chose /eiup = 0.3, as suggested by 
numerical simulations (Walker, Mihos & Hernquist 1996; 
Barnes 1998) which show that the merging galaxy mass 
ratio must be in the range 0.3 to 1.0 to destroy the pro- 
genitor components and produce a spheroidal remnant. The 
value of £*, burst is best constrained by the luminosities of 
bursting galaxies since these are directly related to the star 
formation rate in the burst. To accurately predict the lumi- 
nosities of starbursts a detailed treatment of the absorption 
and reprocessing of light by dust is needed. Such an analysis 
has been carried out by Granato et al. (2000), who found 
that burst = 0.5 provided the best fit to the z — 0, 60/im 
galaxy luminosity function (in which we see light absorbed 
and re-emitted by dust in starbursts) so we will adopt this 
value as default in this work. 

Spheroidals may also form via the gradual accretion of 
smaller galaxies (minor mergers). In the case of a minor 
merger (M1/M2 < /eiUp) we assume no burst of star for- 
mation is involved. Stars from the smaller galaxy are added 
to the bulge of the larger galaxy, which may then become 
dominant, while any residual gas is added to its disk. 

2.2 Model Comparisons based on Number Counts 
and Redshift Distributions 

In this section we are concerned with exploring those param- 
eters central to the taxonomical definition of a morphologi- 
cal spheroidal. The majority of previous studies using semi- 
analytic models to study the properties of galaxies selected 
by their morphological type (Kauffmann 1996; Baugh, Cole 
& Frenk 1996a,b; Kauffmann & Chariot 1998a) have used 
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the bulge-to-total ratio, B/T, usually (as here) measured in 
terms of luminosity. (Specifically, in this work we use the 
dust-extinguished I8i4-band luminosity to define the bulge- 
to-total ratio, B/T;^^^.) Baugh, Cole & Frenk (f996a) also 
considered the time since the most recent major merger 
for each galaxy to be a defining property, placing galaxies 
whose last major merger occurred less than 1 Gyr earlier into 
the Irr/Pec category. Accordingly, spheroidals in the Baugh, 
Cole & Frenk (1996a) scheme jointly satisfy B/T> 0.4 and 
have experienced no major merger in the last 1 Gyr. 

The key issue, central to a successful interpretation of 
the resolved colour data, is that spheroidals can appear to 
be dynamically relaxed while showing evidence of recent star 
formation (e.g. in their internal colours). Instead of adopt- 
ing a fixed relaxation period to define a semi-analytical 
spheroidal, we include a criterion based on the dynamical 
timescale of the remnant. After a major merger, the remnant 
begins to relax from the inside out in a time comparable to 
the local dynamical time (Barnes & Hernquist 1992). For 
each remnant we calculate the dynamical time at the half- 
mass radius. Any object which experienced a major merger 
less than A^roiax dynamical times ago is considered to be an 
Irr/Pec galaxy, irrespective of its B/T. A^reiax is a parameter 
of the model which we expect to be of order unity. For the 
sample of galaxies considered in this work, spheroidal dy- 
namical times are typically a few times 10^ years, i.e. much 
shorter than 1 Gyr. Crucially, for small values of A'^rciax, a 
merger remnant can be classed as relaxed well before the 
distinctive blue colour of recent star formation activity has 
declined. 

We summarise our scheme for defining a sample of mor- 
phological spheroidals in Table ^. Since in this work we are 
primarily interested in the E/SO class only the parameters 
B/Ts and A^roiax need to be constrained. As such, we will 
not atternnt to constrain the values of the other parameters 
in Table |l|. Ahhough Cole et al. (2000) briefiy considered 
the morphological mix of bright galaxies, they did not ex- 
plore in detail the sensitivity of their parameter choices to 
data on morphologically selected faint galaxies. As noted by 
Baugh, Cole & Frenk (1996b) we have some guidance as to 
the value of B/Ts from comparisons of galaxy T- types with 
their bulge-to-total ratio obtained from bulge/disk decom- 
position studies. Such comparisons suggest B/Ts ~ 0.4, but 
with significant uncertainty. Keeping this figure in mind, we 
now explore what constraints can be set on B/Ts and A'rdax 
from the available observational data, specifically morpho- 
logically selected number counts and redshift distributions. 

Number counts are only weakly dependent on the as- 
sumed model parameters. Figure |^ shows the E/SO and total 
number counts in the 78i4 band|]; clearly the baseline model 
(solid lines) matches these data very well. The key question 
is the extent to which the morphologically-classified counts 
depend on our parameters B/Ts and A'reiax. Both the S and 
Irr/Pec number counts (not shown) are relatively insensi- 
tive to the value of A'^^ciax and thus constrain the parame- 
ter B/Ts to be ;^0.44 (in good agreement with the figure 
expected from a comparison to T- types). A value as low 
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Figure 2. Differential number counts in the Igi4 pass-band as 
a function of morphological type. Panels show the counts for all 
galaxies (top) and E/SO galaxies (bottom). Points show a compi- 
lation of observational determinations (as indicated in the upper 
panel) while lines show the model results (solid line — default 
parameters; dashed line — alternative B/Tg; dot-dashed line — 
alternative A^relax)- 



We use Vega-based magnitudes throughout this paper. 
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Table 1. Morphological classification scheme. Galaxies are classi- 
fied on the basis of their bulge-to-total ratio in dust-extinguished 
I-band light, B/Tig^^ (column 1), and the ratio of the time since 
the last burst of star formation to the dynamical time of the 
spheroid, iburst/*dyn (column 2). The morphological class as- 
signed to each region of this two-parameter space is given in col- 
umn 3. 



B/Ti, 



^burstAdyn Morphological class 



0< B/T <B/Tl„ > A^rolax 

B/Ti„< B/T <B/Ts > Af.eiax 

B/Tg<B/T<l >Af.olax 

0< B/T <1 < AT.elax 
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Table 2. Parameters of the relevant morphological classification 
schemes (but which do not affecting the merging and formation 
histories of model galaxies). 

Parameter set 
Parameter Default Alt. B/Ts Alt. A^rdax 



B/Ts 
N 1 



0.44 
3 



0.30 
3 



0.44 
30 



as B/Ts = 0.30 remains consistent with the counts, while 
lower values begin to underpredict the number of spirals. 
The predicted E/SO counts are quite insensitive to A^'rciax 
for A^rciax <10. The dot-dashed line in Fig. ^ shows the re- 
sults for A'reiax ~ 30, which begins to deplete the number of 
E/SO galaxies below that observed. While the counts sug- 
gest a reasonably low value of A^reiax, they do not provide a 
strong constraint. This insensitivity might be perceived as 
an advantage; a detail in the assembly history is irrelevant 
in interpreting an observable. However, as the relaxation 
timescale is central to the rate at which systems become 
recognisable spheroidals, this would be false security, par- 
ticularly in diagnosing systems at high redshift. 

In the hght of the above constraints, in the remainder of 
this paper we will consider three sets of values for the param- 
eters B/Ts and A^'rdax as listed in Table ^. Note that these 
parameters determine only the morphological classification 
scheme and, as such, we would ideally constrain them accu- 
rately from morphologically selected number counts, redshift 
data etc. Later, we will consider altering other parameters 
of the model which directly alter the formation and merging 
histories of model galaxies. It is these formation and merg- 
ing rates which we aim to constrain through the methods 
described in this paper. 

In the upper two panels of Fig. |^ we compare the red- 
shift distribution of model galaxies in magnitude limited 
samples with the recent HST-based redshifts surveys as sum- 
marised by Brinchmann (1999). These encompass the initial 
survey of Schade et al. (1999) enlarged to include a total 
of 54 spheroidals spectroscopically-complete to I8i4 = 22. 
Additional data to a somewhat brighter magnitude limit is 
in principle available from the survey of Im et al. (1999) 
and Im et al. (2001). However, we chose not to include this 
sample as the spheroidals were defined according to sym- 
metry and colour, i.e. differently than in Figure Im et 
al. (2001) discuss the uncertainties of their definition which 



affect the absolute abundances at the ~ ±25% level. As 
expected, the proportion of high redshift spheroidals is sen- 
sitive to A'roiax, with larger values producing fewer examples 
(since more galaxies are placed into the Irr/Pec class). Un- 
fortunately, although the range of redshifts spanned by the 
model galaxies is comparable to that seen in the data, the 
detailed correspondence between theory and observations is 
not particularly good for the total sample (upper panel), a 
point discussed by Baugh, Cole & Frenk (1996a). 

The lower panel of Fig. |^ shows the deeper (I8i4 < 24) 
E/SO sample of Menanteau et al. (1999), which is employed 
extensively in this paper. Barring one discrepant point, the 
model is in very good agreement with this redshift distribu- 
tion. The differences between the model predictions for the 
three parameter sets illustrate again how well these param- 
eters could be constrained by larger datasets. 

While our baseline semi-analytical model, in which 
spheroidals are formed by the mergers of two pre-existing 
galaxies, produces some of the observed trends in this pop- 
ulation, at a quantitative level it disagrees with the present 
redshift data. More relevant at this stage, however, is the 
poor sensitivity of such observables to key physical features 
of the hierarchical model. Even with exquisite data at these 
magnitude limits, we can expect it will be hard to conclu- 
sively determine whether the hierarchical model is correct 
given the freedoms in the morphological definitions intro- 
duced in Table 1. Therefore, we propose to use these ob- 
servational data to constrain the morphological classifica- 
tion parameters of Table ^ and then use other data, such as 
that described in this paper, to constrain the formation and 
merger rates of spheroidals. 



2.3 Quantifying Recent Activity Using 
Colour-Based Properties 

In the present section we demonstrate that the colour in- 
homogeneities discovered for HDF spheroidals (Menanteau, 
Abraham & Ellis 2001) offer a much more direct probe of 
the hierarchical hypothesis tha n th e number count and red- 
shift distributions discussed in §2.2, Not only do the internal 



colour variations seen in the HDF offer an important way 
to separate models based on isolated passive evolution fol- 
lowing monolithic collapse from those invoking continuous 
hierarchical assembly, but they also provide the first oppor- 
tunity for a quantitative measure of the assembly rate in 
the latter case. Importantly, Menanteau, Abraham & Ellis 
(2001) find that around one third of spheroidals in their sam- 
ple show marked inhomogeneities in their internal colours, 
and that in the majority of these cases the inhomogeneity is 
characterised by a blue core in the galaxy colour profile. The 
observations therefore require that blue light be more cen- 
trally concentrated than red light. Here we will assume that 
these blue cores arise due to recent bursts of star formation. 

For present purposes, where we are concerned only with 
V and I photometry, the major merger event which forms 
each elliptical in the model can be characterised by two phys- 
ical quantities: (i) the fraction of stars formed in the associ- 
ated burst relative to the final stellar mass of the spheroidal 
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Figure 3. (cont.) The three panels show faint galaxy redshift 
distributions, with observational data indicated by filled circles 
and model predictions shown as histograms. In the upper panel 
we show the redshift distribution of galaxies of all morphological 
types brighter than I814 = 22 from Brinchmann (1999), while the 
middle panel shows the redshift distribution of spheroidal galaxies 
from the same sample. The lower panel shows the redshift distri- 
bution of Igi4 < 24 spheroidals from Menanteau et al. (1999). In 
the lower two panels, model results are shown for standard (solid 
histogram), alternative B/Ts (dashed histogram) and alternative 
A^reiax (dot-dashed histogram) parameters. 



remnant, /burst, and (ii) the time since the merger, tburstQ- 
The first determines how the colour of the remnant evolves 
with time, and which can in principle be constrained by the 
resolved HDF data, and the second determines at what point 
in that evolution the galaxy is observed. Ideally we would 
like to measure these two parameters directly from the ob- 
servational data as this would allow us to reconstruct the 
rate of spheroid formation as a function of redshift. However, 
with the present photometric data and relatively small sam- 
ple size such an approach is not warranted (spectroscopic 
measurements will provide far better constraints on these 
parameters). Instead we focus on the much simpler objec- 
tive of assessing the compatibility of hierarchical formation 
with the present data. However, we will give predictions for 
the variation of these key parameters with redshift which 
may then be tested by forthcoming observations. 

In Fig. ^ we show how the two physical quantities vary 
with redshift for a magnitude- limited sample of model galax- 
ies in a ACDM cosmology. In this and subsequent figures, we 
will consider spheroidal galaxies brighter than /gi4 = 24 in 
order to match the HDF sample used by Menanteau, Abra- 
ham & Ellis (2001). Unless otherwise stated all model re- 
sults refer to our default parameters (refer to Table ^) . 

There is little variation in /burst with redshift, but a 
significant dispersion at fixed z. The median burst typically 
involves 10-20% of the total spheroidal mass. The lack of any 
trend with redshift is the result of two counteracting effects. 
Firstly, /burst shows a strong correlation with the mass of 
the spheroidal galaxy, such that lower mass galaxies tend 
to have higher values of /burst- This trend arises because 
more massive ellipticals tend to form from more massive 
disk progenitors Kauffmann & Chariot (1998a) and, in the 
model of Cole et al. (2000), more massive disks tend to have 
smaller gas fractions (at least if we consider only central 
galaxies, which are always involved in the formation of an 
elliptical), a consequence of the variation of star formation 
rate with disk mass. 

The apparent magnitude limit of the sample prevents 
us from seeing the least massive (highest /burst) galaxies at 
high redshifts. However, at high redshift there are also fewer 
of the most massive (lowest /burst) gala^cies (since these have 
not had time to form). 



Figure 3. ^ The metallicities of pre-existing and newly-formed stars will 

also have some effect on the colours, an effect correctly accounted 
for in the model, but is of secondary importance to the present 
discussion. 
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Figure 4. The variation with redshift of the two parameters defining a merger event, /burst a-nd ^burst (left and right-hand panels 
respectively). Lines show the median of the model relation with error bars indicating the 10% and 90% intervals. Solid lines show results 
for our default parameters, dashed and dot-dashed for the alternative B/Tg and iV^eiax respectively. For clarity, error bars are shown 
only on the solid lines. The 10% and 90% intervals for the other models are comparable to the default model with the exception that 
the 10% intervals of fburst for the alternative A'^^clax model occur at significantly larger tburst- 



There is a strong trend of tburst with redshift however. 
This quantity is bounded by the age of the universe at given 
z (since no galaxy can have experienced a burst which began 
earlier than t=0). The dotted line in Fig. ^ shows the age of 
the universe as a function of redshift. The evolution in the 
median value of fburst is much more rapid than this. This 
is a consequence of higher rates of merging at high redshift 
combined with the fact that, at high redshift, a galaxy is seen 
only if it is intrinsically luminous and so we preferentially 
select galaxies which have recently experienced a burst of 
star formation, these being more luminous than quiescent 
examples (particularly in view of the k-correction at z m 
1.5). 

In conclusion the model predicts that for the magnitude 
limited sample considered here, a strong correlation between 
iburst and redshift should exist, such that higher redshift 
spheroidals have typically experienced a burst more recently. 

We now describe how we measure the degree of recent 
activity in a spheroidal galaxy from the available photo- 
metric data (both real observational data and mock images 
of model galaxies). Menanteau, Abraham & Ellis (2001) 
proposed a statistic 5(V — I) which measures the degree of 
inhomogeneity in the internal colours defined as: 



S(y -l) = 2N 



(1) 



where Xi is the colour of the i**^ pixel, SNR(a::i) is the 
signal to noise ratio, S{xi) is a step function which includes 
only those pixels with SNR(a;i) > 1.3 in the sum and N is 
an arbitrary constant which was taken to be A'^ = 3. The 
sum is taken over all pixels brighter than Veoe = 25.5mag 



arcsec"^. The value of (5(V — I) obtained for any single galaxy 
will depend on the point spread function (PSF) of the ob- 
serving instrument, since substantial smoothing will tend to 
homogenise the internal colours. Although wavelength de- 
pendent PSFs, such as occur in the WFPC2, can create ar- 
tificial colour inhomogeneities. Menanteau, Abraham & Ellis 
(2001) found this to be a negUgible source of inhomogene- 
ity for the HDF spheroidals except for a few galaxies with 
highly concentrated surface brightness profiles. 

As an alternative diagnostic we consider the differ- 
ence in integrated colour compared to that expected at the 
galaxy redshift for an old, passively evolving population, 
(V — I)passivG(z). The latter case is defined here as a sys- 
tem which formed at = 3 with an exponentially decaying 
star formation rate with e-folding time of 1 Gyr and with 
Solar metallicity. In Fig. ^ we show the integrated V — I 
colours of observed and model spheroidal galaxies as a func- 
tion of redshift and compare these to (V — I) passive (2:). We 
plot data from the HDF sample of Menanteau, Abraham & 
Ellis (2001) and also from the brighter (Iab < 22) sample 
of Schade et al. (1999) to demonstrate that both datasets 
display similar trends. Note that while both model and data 
broadly trace out the shape of the passive curve, both reveal 
many galaxies that are significantly bluer. This is particu- 
larly noticeable for the model at z ^1 where the median 
colour is ~ 0.3 magnitudes bluer than the passive predic- 
tion. 

Interestingly, the HDF data also contains a few very red 
objects (i.e. redder than the passively evolving track shown 
in Fig. bl and in fact redder than passively evolving tracks 
for galaxies formed even at higher zf)- The most anoma- 
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Figure 5. Integrated V— I colour as a function of redshift. Cir- 
cles show data from the HDF North and South with solid points 
indicating spectroscopic redshifts and open points photometric 
redshifts. Triangles show brighter data from the survey of Schade 
et al. (1999) (with morphologies classified visually). The median 
of the model distribution is shown by the solid line, with error 
bars indicating 10% and 90% intervals of the distribution. The 
dashed line indicates the colour of a galaxy formed at z = 3 with 
an exponentially decaying star formation rate with e-folding time 
of 1 Gyr and with Solar metallicity. 



lous sources are ones with photometric redshifts above unity. 
The red colours could arise as a consequence of an overes- 
timated photometric redshift although two of these galaxies 
would be extremely red whatever redshift they were placed 
at. Certainly, Fig. ^ suggests that our present model does 
not produce spheroidal galaxies as red as those observed in 
the HDF. 

We define the colour excess: 



(V-I)c 



= (V-I)-(V-I)paBsivc(2). 



(2) 



While this simple measure has the desirable property of be- 
ing largely independent of the PSF it makes no use of the 
resolved internal structure, and furthermore is sensitive to 
errors in the redshift of the galaxy (which is particularly 
important for the galaxies with only photometric redshifts) . 
Nevertheless, it is interesting to compare this measure with 
the 5iV - I) index. 

To determine the sensitivity of both statistics for the 
model galaxies we construct a mock image of each galaxy 
(including noise which makes a small contribution to 5(V — 
I)), smoothed with a Gaussian PSF of FWHM 0.126 and 
0.143 arcseconds in the V and I-bands respectively to mimic 
the PSF of the WFPC-2 (Casertano et al. 2000). Both 
statistics described above are then measured as for a real 
galaxy. For the purposes of constructing mock images each 
galaxy is assumed to consist of three components, a disk, a 
spheroidal remnant and a recent burst. Although a galaxy 



may have experienced more than one burst in its history we 
treat only the stars formed in the most recent burst as a 
separate component, with stars which were formed in ear- 
lier bursts assumed to be uniformly mixed in the spheroidal 
component. Our results are insensitive to this assumption, 
since placing the light from stars formed in all the bursts in 
a separate component does not make any significant differ- 
ence to the statistics we consider.) Disks are modelled with 
an exponential surface brightness profile and random incli- 
nation. Hernquist (1992) has shown that mergers between 
purely stellar systems, specifically bulge-less galaxies, result 
in a remnant which has a surface brightness profile close to 
an r^'''*-law, but with a core of constant surface brightness. A 
simple parameterisation which displays qualitatively similar 
behaviour is 



E(r-) = Eoexp(-7.67 {[r/r,] 



l2 , 2 
+ C 



1/8 



(3) 



where r is the radius, re the effective radius and c is a core 
radius in units of re. We will refer to this form as the "Hern- 
quist profile" hereafter. This expression is not proposed as 
an accurate fit to the Hernquist (1992) results, but has 
the desirable properties of being analytically integrable and 
able to produce remnants with any desired core size. We 
distribute the mass of the pre-existing stars with this sur- 
face density profile. Gas in mergers is generally driven to 
the centre of the remnant where it presumably undergoes 
a burst of star formation (Barnes & Hernquist 1996). The 
stars formed in the recent burst are therefore distributed in 
such a way that the total surface mass density (i.e. the sum 
of burst and pre-existing components) follows an r^'^^-law 
(i.e. as eqn. ^ with c = 0). The total mass in the Hernquist 
profile is given by M(c) = 27rrE(r)dr. Requiring that 
M(c)/M(0) = Mpi.o/(Mburst + Mpro), where Mburst is the 
mass of the stars formed in the burst and Mpre is the pre- 
existing stellar mass, allows us to solve for c. The value of 
Tc is determined from the half-mass radius of the remnant 
as calculated by the galaxy formation model (Cole et al. 
2000). The light of each component is then distributed as 
the mass assuming a constant mass-to-light ratio. With this 
construction the spheroid would eventually come to resemble 
an r'^''*-law at late times when the stars of both components 
were old (ignoring any differences in metallicity between the 
two populations). 

In cases where the recent burst accounts for a large frac- 
tion of the total mass of the spheroid the core radius, c, can 
become very large. No simulations of mergers between such 
gas-rich disks have been carried out (but see Somerville, Pri- 
mack & Faber 2001), so the above form for the pre-existing 
and burst star populations may not be relevant in such cases. 
We therefore consider an alternative case in which both pre- 
existing and burst stars are distributed as r'^''*-law profiles, 
but with the burst stars having a much smaller effective ra- 
dius (one tenth of the spheroid radius) . 
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3 ANALYSING RECENT ACTIVITY IN HDF 
ELLIPTICALS 

We now compare the distributions of model and observed 
HD F sp heroidal galaxies using the two statistics described 



in §2.3 



The lower left-hand panel of Fig. ^ shows differential 
distributions of each statistic. For (V — I) excess the model 
produces a comparable distribution to the data for blue 
objects ((V — I)exeess ^0), but caunot account for the red 
objects ((V — I)excess >0) as was expected on the basis of 
Fig. ^. This quantity is independent of the surface bright- 
ness profiles we ascribe to model galaxies, but does depend 
on the values of the parameters B/Ts and A'reiax- In the 
lower right-hand panel of Fig. ^ we demonstrate the effects 
of using the alternative parameter values given in Table ^ 
Using B/Ts = 0.30 instead of 0.44 produces substantially 
more galaxies just bluewards of (V — I)excess = (this al- 
ternative morphological definition tends to include several 
galaxies which experienced their last burst much longer than 
IGyr ago and whose only star formation is occurring in their 
disks), without affecting the rest of the distribution. Increas- 
ing A'^reiax from 3 to 30 significantly reduces the number of 
very blue model galaxies. With the larger value of A'reiax 
these extremely blue galaxies, which have typically experi- 
enced a major merger in the very recent past, are rejected 
from the E/SO class. With the present data we can infer 
A'reiax <10. Notc that neither of these parameter alterations 
helps reconcile the observed and predicted distributions of 
red galaxies. 

For ^(V — I) we show two model lines in the upper left- 
hand panel, the solid one using the Hernquist profile and the 
dashed one using an r^'^^-law for both burst and pre-existing 
stellar components. The choice of profile does not drastically 
alter the width or shape of the distribution. Using r^''*-law 
profiles produces somewhat more galaxies at low values of 
S(V — I) and does not produce the tail to high values which 
occurs with the Hernquist profile. We remind the reader that 
while the Hernquist profile has some theoretical justification, 
the r^'^^-law profile has none. The purpose of considering this 
profile is to demonstrate that the model predictions are not 
overly sensitive to this choice. 

In the upper left-hand panel we again see that the model 
distribution is in reasonable agreement with the data for 
high values of 5(V ~ I), but appears to predict too few 
galaxies with low 5(V — I). In the upper right-hand panel 
we consider the effects of changing model parameters on the 
predictions for 5(V - I). Changing B/Ts from 0.44 to 0.30 
produces more model galaxies at low (5(V — I), thereby im- 
proving the agreement between model and data (although 
some differences remain), while increasing A'^reiax to 30 de- 
pletes the number of galaxies at all 5{Y — I), but partic- 
ularly so at high S(V — I) since these galaxies are those 
which have experienced a major merger in the very recent 
past. However, even with this large value of A'reiax (which is 
marginally ruled out by the morphologically selected counts 
of Fig. ^) the model still provides a reasonable match to the 
high (5(V - I) data. 

Finally, in Fig. ^, we compare model and data as a func- 
tion of redshift. For the model the median relation is shown 



by the line, with error bars showing the 10% and 90% in- 
tervals at each redshift. Error bars are typically shown on 
only one line per panel to aid clarity but are of comparable 
size for the other lines. Data are shown by circles; filled and 
open to indicate spectroscopic and photometric redshifts re- 
spectively. As shown in the upper panels, both model and 
data show an increasing median 5(V — I) with redshift out 
to z « 1. The top left-hand panel clearly demonstrates that 
the model predictions for 5(V — I) at high redshifts depend 
crucially on the surface brightness profile adopted for the 
post-merger spheroidal, with the Hernquist profile produc- 
ing many more high S(V — I) galaxies at high redshift than 
the r^^^ profile (for which the distribution of 5(V — I) is in 
rather good agreement with the data at all redshifts) The 
exact choice of model parameters has little effect on the 
prediction (although increasing A'reiax to 30 does reduce the 
90% intervals). 

In (V — I)exeess (lower panels) the model spans a similar 
range to the data as a function of redshift, the rapid increase 
in blue excess objects between z — and z = 0.5 is apparent 
in both model and data for example. Interestingly, the data 
appears almost bimodal at these redshifts, with one group of 
galaxies lying along (V — I)excess ~ (and so are presumably 
old, passive spheroidals) and another group lying along a 
locus of (V — I)exeess which dccrcascs with redshift. Dots 
show a sample of model galaxies for comparison. Although a 
clear branch of model points can be seen running along (V — 
I)excess ~ 0, uo well-dcfiued branch is seen in the vicinity of 
the second group of data points, although the model points 
certainly span the range of blue (V — I) excess values seen 
in the data. Given the small number of data points (and 
the fact that most of those in the second branch have only 
photometric redshifts) no strong conclusions can be drawn 
as yet. Of course, the model does not contain the very red 
galaxies seen in the data, as expected from Fig. |^. 

Since we can easily switch off the effects of the PSF on 
model galaxies we have checked that the wavelength depen- 
dence of the PSF does not significantly affect the distribu- 
tion of S(V — I) for model galaxies. Without the PSF the 
distribution of 5(V — I) for model galaxies is significantly 
higher in the lowest bin of (5(V — I) shown in Fig. ^ (the 
wavelength dependence of the PSF increases 5(V — I) for a 
few of these galaxies, mostly at low redshifts, giving them 
5{V-l) « 0.1) and the taU extends to slightly higher 5(V-I) 
(here it is the homogenising effects of the PSF which are im- 
portant, particularly for high redshift galaxies). 

We have also checked the effect of the disk compo- 
nent of the model spheroidal galaxies. Recall that our E/SO 
class contains galaxies with at least 44% of their light in a 
spheroidal component, allowing up to 56% to be in a disk 
(although many have much less than this of course). The 
disk light can contribute to S(V — I) and (V — I)excess since 
it may be of a different colour to the spheroidal light and is 
also distributed differently. The broad features (e.g. width) 
of the 5(V— I) and (V — I)execss distributions remain the same 
whether we include the disk components of E/SO galaxies or 
not. 

Without the disk component we find that the tail of 
the S(V — I) distribution extends to slightly higher values 
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Figure 6. Differential distributions of both statistics discussed in §2.3 for E/SO galaxies brighter than I814 = 24. Circles show observa- 
tional data from the combined HDF North & South. Histograms show model results. In the left-hand column solid histograms are for the 
Hernquist profile, while dashed histograms use an r^/'*-law profile. In the right-hand column solid lines are for the default parameters, 
dashed for for alternative B/Tg and dotted for alternative A^rclax- 



(in galaxies with a recent burst adding a disk component 
can actually reduce (5(V — I) since the disk may be of sim- 
ilar colour to the burst component), but the effect is small 
relative to the current theoretical and observational uncer- 
tainties. In (V — I)excess leaving out the disk components re- 
duces the number of galaxies at (V — I)oxccss ~ —0.4 to —0.1, 
since model disks often have (V — I)cxccss in this range. In 
conclusion, including the disk component of our spheroidal 
galaxies has some effect on the predicted distributions, but, 
with the present dataset, our basic conclusions are insensi- 
tive to the presence or absence of a disk component in model 
galaxies. 

So far we have considered most of the parameters of 
the model to be fixed at the values chosen by Cole et al. 
(2000), the exceptions being B/Ts and A'rciax which we have 
constrained from morphologically selected number counts 



and redshift distributions. By adopting this philosophy the 
model results are true predictions, since the model param- 
eters were fixed by considerations entirely unrelated to star 
formation activity in spheroidal galaxies. 

Finally, we now consider the effect of altering key model 
parameters which control the formation and merging histo- 
ries of model galaxies. We will not explore this graphically 
but it is interesting to examine the consequences qualita- 
tively for two reasons. Firstly, this will aid our understand- 
ing of what physical processes are crucial for constraining 
the statistics we have considered. Secondly, and most im- 
portantly, it will test the ability of the type of compari- 
son between theory and observation developed in this pa- 
per to place useful constraints on the formation rates of 
spheroidal galaxies. An important caveat must be borne in 
mind — namely that after we alter these parameters it is 
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Figure 7. The variation of our two statistics with redshift. Top and bottom rows show 5(V — I)and (V — I)cxccss respectively. In each 
case circles show HDF galaxies (filled points are spectroscopic redshifts, open points are photometric redshifts). Lines show the median of 
the model distribution with error bars indicating the 10% and 90% intervals. In the left-hand column the solid line shows results for the 
Hernquist profile (the dots show a sample of model galaxies to give a more detailed indication of the distribution of model results), while 
the dashed line shows the result for the r^/''-law profile. In the right-hand column the solid lines show results for our default parameters, 
while dashed and dotted show results for the alternative B/Ts and alternative A'^rolax respectively. 



unlikely that the model will be in as good agreement with 
the morphologically selected counts and redshift distribu- 
tions of Figs, y & H as with the standard set of parameters. 
Such models are not necessarily physically realistic, and we 
explore them only to demonstrate the sensitivity of 5(V — I) 
and (V — I)cxcoss to the formation histories of spheroidals. 

Three model parameters are crucial to the formation 
of spheroidals in the model of Cole et al. (2000). The fre- 
quency of mergers is controlled by the parameter /df which 
determines the strength of the dynamical friction force felt 
by galaxies. The division between minor and major mergers 
is determined by the parameter /oiup and the timescale of 
the star burst associated with a major merger depends upon 



£*, burst- (See and Cole et al. (2000) for further details.) 
Increasing /df enlarges the merger timescale and so reduces 
the frequency of the major mergers. Increasing this param- 
eter from 1 to 2 or 4 substantially reduces the number of 
spheroidal galaxies. The fractional reduction is greatest for 
high (5(V — I) (highly negative (V — I) excess) galaxies. This 
arises since the reduced rate of major mergers means that 
a given spheroid tends to have experienced its last burst of 
star formation longer ago, which results in the characteristic 
colours of the burst stars fading away before the galaxy is 
observed. 

Increasing /eiup to 0.5 or 0.75 also has a large effect on 
the distributions of 5(V — I) and (V — I)excess. Firstly, there 
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is an overall reduction in the number of spheroidals, since 
major mergers become much rarer. There is also some diflFer- 
ential effect, with high (5(V — I) (equivalently, highly negative 
(V — I)excess) galaxies being depleted preferentially, again be- 
cause of the greater time since the most recent merger. A 
secondary effect of the reduced merger rate is that, partic- 
ularly at low redshift, galaxies in major mergers tend to 
be less gas rich when /dup is increased, as they have had 
much more time to turn gas into stars quiescently. This also 
tends to reduce the number of high (5(V — I) (highly negative 
(V - I)oxcoss) galaxies. 

The parameter e*, burst has little effect on the distribu- 
tions of (5(V — I) and (V — I)cxccss since it does not change 
the frequency of bursts nor the mass of stars formed within 
them. It does however alter the sensitivity of the results to 
the value of A^reiax- For smaller e*, burst we find that increas- 
ing A^reiax produces a greater reduction in the numbers of 
high 5(V — I) and highly negative (V — I)excess galaxies. This 
occurs because, with the slower star formation rates implied 
by a low e*, burst, bursts tend to redden more rapidly and so 
become indistinguishable from the rest of their galaxy. With 
A'reiax = 3 many model spheroidals still have recognisable 
blue colours to their cores, but for A^reiax = 30 most have 
already reddened sufficiently to deplete the highly negative 
regions of the (V — I)cxccss distribution. 

The Cole et al. (2000) model assumes that in the case 
of a minor merger (M1/M2 < /ciup) stars from the smaller 
galaxy are added to the bulge of the larger galaxy, thereby 
providing a mechanism by which bulges can form without 
any associated burst of star formation. An alternative would 
be to assume the stars are added to the disk of the larger 
galaxy, in which case spheroids can only ever form through 
a major merger. Cole et al. (2000) found that the choice of 
where to place stars from minor mergers made little differ- 
ence to the properties of model galaxies which they consid- 
ered. We find a similar result, with no significant difference 
being made to the (5(V — I) and (V — I)oxccss distributions or 
their redshift dependence. 



4 DISCUSSION 

We have assessed the implications of the sample of elliptical 
galaxies in the HDFs with inhomogeneous internal colours 
(Menanteau, Abraham & Ellis 2001) for the picture of con- 
tinuous formation of ellipticals expected in hierarchical mod- 
els of structure formation. 

We have shown that the galaxy formation model of Cole 
et al. (2000) can produce good agreement with data on mor- 
phologically selected number counts, and produces many of 
the trends observed in morphologically selected redshift dis- 
tributions, although the detailed agreement in these quanti- 
ties is less good. In defining the morphology of model galax- 
ies we considered the time taken for a merger event to relax 
to an elliptical galaxy configuration, quantifying this cri- 
terion in terms of a multiple, A'roiax, of the elliptical's dy- 
namical time. Of crucial importance to the present work is 
that this definition allows a galaxy to be classed as a re- 
laxed elliptical before the distinctive blue colour of recent 



star formation has faded. With the present level of agree- 
ment between theory and observations it is difficult to set 
strong constraints on N^^ia,^, but A^'rdax ;$10 is favoured, as 
would be expected on theoretical grounds (Barnes & Hern- 
quist 1992). Accurate determinations of the redshift distri- 
bution of elliptical galaxies will provide the strongest con- 
straint on this parameter. 

We compare the observational data with model predic- 
tions through two statistics, the inhomogeneity index origi- 
nally defined by Menanteau, Abraham & Ellis (2001) and 
the difference between the observed (V-I) colour and that 
of a passively evolving old stellar population at the same 
redshift (which is simpler to compute from the model, but 
makes no use of the spatially resolved information in the 
HDF images and which is sensitive to errors in the photo- 
metric redshifts of HDF galaxies). 

The theoretical predictions suffer from some uncertain- 
ties, due to the rate at which ellipticals are expected to form 
in a given cosmological model (as demonstrated in Fig. 
the uncertain spatial distribution of light from recent star 
formation for which we have only limited theoretical guid- 
ance as yet and the parameters used to specify what consti- 
tutes a spheroidal galaxy in the model. Nevertheless, we have 
shown that our basic conclusions are insensitive to the latter 
two uncertainties (the first is exactly what we are attempt- 
ing to measure from the data so necessarily has a strong 
impact on the model predictions), allowing us to make a 
useful comparison to the observational data. 

Comparing our two statistics for model and observed 
galaxy samples we find reasonable agreement in the distri- 
butions of both (5(V — I) and (V — I)cxccss with two no- 
table exceptions. Firstly, the current model produces too 
few spheroidals with low 5(V — I) and, secondly, the model 
does not produce the population of (V — I)oxcoss > galaxies 
seen in the data. We discuss these two points further below. 
However, most importantly for the present work is the fact 
that the model does produce the correct number of high 
5{y — I) and blue (V — I)cxccss galaxies. The implication of 
this result is that a model in which elliptical galaxies form 
hierarchically in a ACDM universe, and which was specif- 
ically constrained to match a set of local data (Cole et al. 
2000), is quantitatively consistent with the continued star 
formation implied by the HDF data of Menanteau, Abra- 
ham & Ellis (2001) with no further adjustment of the free 
parameters of the model. The basic trends of 5(V — I) and 
(V — I)cxccss with redshift are also reproduced by the model. 
An expanded data set will allow these trends to be tested 
quantitatively. 

However, there is evidence that the model produces 
somewhat too few relaxed, homogeneous spheroidals from 
the 5{y — I) statistic. Given the small sample size and the 
theoretical uncertainties it is premature to assign much sig- 
nificance to this apparent discrepancy. Perhaps the most no- 
ticeable discrepancy between theory and observations is the 
lack of extremely red objects (relative to a passively evolv- 
ing spheroidal formed at high redshift) in the model sample. 
The majority of these objects have only photometric red- 
shifts and so it remains possible that their apparently very 
red colours for their redshift are simply a consequence of 
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an incorrect redshift determination. Both of these points 
should be addressed in more detail, both using larger ob- 
servational samples, and by improvements in the theoretical 
calculations. 

We have briefly considered the sensitivity of the 5(V — 
I) and (V — I) excess statistics to the formation rate of 
spheroidals, by varying parameters of the model which di- 
rectly influence this formation rate. We typically find that 
lower formation rates tend to preferentially deplete the high 
5{V — 1) and highly negative (V — l)cxccss regions of the 
distributions. This is due mainly to the fact that, for a 
lower spheroidal formation rate, a galaxy's last major merger 
tends to have occured longer ago and so the characteristic 
colours of the recent burst have faded away by the time it is 
observed. Thefore, analysis of large samples of spheroidals 
using these techniques should allow us to directly constrain 
their formation rate. 

Spectroscopic observations of these faint populations 
of ellipticals should allow estimates of their star forma- 
tion rates and the ages of their stellar populations. These 
quantities, for which we have presented predictions, can be 
compared to the theoretical predictions much more directly, 
and so should provide much stronger constraints on the 
formation of ellipticals. With the imminent launch of the 
Advanced Camera for Surveys on Hubble Space Telescope, 
there is a strong case for assembling a large sample of field 
spheroidals with redshifts and reliably-sampled internal col- 
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